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Abstract 

The radiative corrections involved in the precise determination of the electroweak mixing 
angle measured at the Z° peak are reviewed in detail, with particular emphasis on the 
new calculation of two-loop heavy top effects. This example serves as a brief pedagogical 
introduction to the problems and techniques of higher order electroweak calculations. 
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The experimental verification of the Standard Model (SM) has reached in the last 
decade a high degree of sophistication. On the experimental side we have now a large 
amount of very precise data collected at high energy, and several observables are known 
at the permille level. On the theoretical side, a major effort has been undertaken in 
order to match the experimental accuracy. The one-loop corrections to all electroweak 
observables are by now very well established H , and two and three-loop effects have been 
investigated in many cases. Because of the number of different mass scales involved, multi 
loop calculations in the SM are highly non-trivial, but in the last few years a surge of 
activity in this field has made the investigation of some higher order effects possible . 

Although the assessment of the theoretical error is a very subjective matter, we can 
generally distinguish between two different kind of uncertainties: parametric, i.e. due to 
the uncertainty of the input parameters, and intrinsic, that is inherent to the truncation 
of the perturbative expansion, or to the appearance of non-perturbative effects. An 
important example of parametric uncertainty is the one due to the hadronic contribution 
to the running of the electro-magnetic coupling from low energy to the Z° resonance. As 
long distance dynamics is involved in this case, one has to resort to the experimental data 
for e + e _ — > hadrons. Through the use of dispersion relations, a(M z ) is then determined 
with an accuracy of about 7 x |J. This sizable uncertainty is the main source of 
theoretical error for the Z°-peak observables (except for low-angle Bhabha scattering). 
On the other hand, scheme dependence, scale dependence, and explicit evaluation of 
higher order corrections can all give information on the intrinsic error, hence the crucial 
importance of higher order corrections for the tests of the SM. 

In this lecture I will use a relatively simple example to illustrate some of the problems 
and techniques involved in the evaluation of the higher order effects relevant for the 
precision tests of the SM. The case of the effective Weinberg angle, sin 2 9^j, measured 
at LEP and SLC, has both the relevance and the simplicity to make a detailed and 
self-contained discussion possible. This observable is now known with excellent accuracy 
(the present world average is sin fljjfj = 0.23165 ±0.00024 |§), and its sensitivity to the 
Higgs mass is very high in comparison to other quantities, giving it a predominant role 
in the present analyses. 

After introducing to the basics of the renormalization of the SM, I will review the 
one-loop electroweak corrections to sin 2 d^ff ■ I will then discuss some potentially large 
higher order effects: QCD corrections to hadronic loops, large mass two-loop effects (top 
quark and heavy Higgs case), etc. In the case of the two-loop heavy top corrections, I 
will describe in detail a new calculation of the 0(a 2 mf/M^) contributions, and illustrate 
explicitly the heavy mass expansion method. 

In the SM Lagrangian, besides the masses of the fermions and of the Higgs boson, 
and the CKM mixing parameters, there are only three independent parameters, which 
we can choose to be, for example, 

9 9' v, or g M w M z , (1) 

where g and g' are the SU{2) L and U(l)y couplings, and v is the vacuum expectation 
value of the Higgs field. On the other hand, there are presently three very obvious 
experimental inputs that we can use to determine these parameters: a, the fine structure 
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constant, G^, the Fermi constant, and M z ; they are respectively known within 4 x 10~ 9 , 
8 x 10~ 6 , 2 x 10~ 5 . What effectively enters the physics of the weak interactions is generally 
a(M z ), the running electro- magnetic coupling at the Z° scale. For the reasons mentioned 
above, this is known only slightly better than the permille level. 

A number of natural relations link the couplings and the masses of the gauge bosons 
occurring in the bare SM Lagrangian. Although initially defined in terms of the gauge 
couplings, after spontaneous symmetry breaking the electroweak mixing angle (P w relates 
both bare masses and couplings among themselves: 

tan 91 = 4 e = g sin 6° w , M° w = cos 6° W M Z . (2) 

At any order in perturbation theory we can relate the parameters of the Lagrangian to 
the measured inputs, calculate them, and then make predictions for any observable. For 
instance, using s = sm9 w and c = cos9 w , most relevant among such relations are 

2 7ra 2 7la 1 fry \ 



M 2 w = c I M I z - M l w = c R M l z p R . (3b) 

Here the l.h.s. corresponds to a tree level description, while on the r.h.s. I have shown 
that, in a given renormalization scheme "R", the quantum effects can be incorporated 
through the radiative corrections Ar* and p R . They are functions of s R , M w , m t , M H , 
etc., and clearly depend on the renormalization scheme. In the following I will adopt 
MS renormalized couplings at the Z° scale, so that the mixing angle is defined by 



- sin w (M z ) m = — (4) 

0£2{M z )jfis 



(a and ai are the MS U(l) e . m . and SU(2) running couplings), while the vector boson 
masses are defined as the physical masses, as in Eq. (|3|) ||. This choice has some 
advantages which will be clear later. In this frameworkF], Ap = 1 — /T 1 and Ar w are 
the radiative corrections entering the l.h.s. of Eq. (|3|). By solving these two equations 
simultaneously, the mass of the W boson and the MS mixing angle can be determined. 

At LEP and SLC sin 2 O^Jl is measured from the on-resonance asymmetries. Leaving 
aside the photon mediated amplitudes and the QED radiative corrections, which form a 
finite and gauge invariant set and are routinely subtracted by the experimental groups, 
the left-right (LR) and forward-backward asymmetries depend only on the ratio of the 
vector and axial-vector couplings of the Z° to the leptons. In the simple case of the LR 
asymmetry measured at SLC, 



■ _ g e A r 



1 Alternative approaches are reviewed in Q; the on-shell scheme was originally proposed in [Q]. 
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where gy and g A are the vector and axial-vector couplings of an on-shell Z° to the 
leptons. 

At the tree level, the amplitude for the decay of a Z° boson on mass-shell is given by: 




where e M is the polarization vector of the Z , and J% is the fermionic current coupled to 
the Z°. The latter can also be written as 

J£ = \ J z- s " J " . ( 6 ) 

with J!* = Q q qj^q, and J3 = I3 q-y^a^q the U(l) and SU{2) L currents respectively 
(If = ±1 and a T = ~(1 =f 75)). We then define 1 — 4sin 2 O^ff = g v /g e A , and Alr directly 
measures sin 2 1 ^. 

At the one-loop level much of this simplicity is retained. We need to consider only 
the electroweak effects; the two topologies that contribute to the Z-decay amplitude at 
this order are depicted in FigJI], where (a) summarizes all vertex corrections and the 




Figure 1: One-loop Feynman diagrams contributing to Z° decay amplitude. 

wave function renormalization of the external fermions, and (b) represents a self-energy 
insertion, that either mixes the Z with the photon, or contributes to the wave function 
renormalization of the Z°, in which case it is multiplied by a factor 1/2. Just because 
of their tensorial structure, the sum of these contributions and of the tree level one can 
be written in terms of bare quantities as 

-i^A(££\ l^-slBJ^\0), (7) 
c z 

where A and B are two divergent quantities, regulated in n dimensions. At this point 
we perform the renormalization of the couplings by a MS subtraction at the Z° scale. 
In Eq. (0) we replace so by s(fi) — 5s, and eo = go so by e(/i) — 5e, and set the scale 
/i = M z . Notice that 5s and 5e are pure poles in e = (4 — n)/2. Expanding everything 
up to 0(a) we get rid of the divergences. The new amplitude is finite: 

-if Vp{^\\j3"-s 2 k(M 2 z )Jt;\0), (8) 
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and we can identify 

sin 2 9 l e f} = Re k{M 2 z ) s 2 (M z ). (9) 



In order to actually compute sin 2 , we therefore need the form factor k(M 2 ), 
which depends only on the component of the diagrams in Fig.|l]. This is given by the 
mixing amplitude, — igs (££\ Jif|0) A^MS^/M 2 ., and by the part of the vertex amplitude 
proportional to J^, —ig/c (££\Ji*\0) s 2 V 1 . Hence 



i(Mj) = 1 _|A^)_ KT _w (10) 



From the ultraviolet divergent part (see for ex. the second paper in Ref. 0) of the vector 
boson mass counterterms we calculate the counterterm 5s 2 = —5c 2 = —[5(Ml/M 2 )]uv, 

^ 9* (K . 11 S A 1_ (n) 



s 2 16tt 2 V 6 ' 3 " J e 
In the t'Hooft-Feynman gauge the divergent part of the mixing self-energy is 

A, tZ (M 2 7 ) g 2 s 1 \7 17 „ Ml - ^. . 

6 + y s +2 #~ 2a +0(e) 



(12) 



M| 16tt 2 c e 

Similarly, the divergence of V~ is 2 g 2 c 2 /(167r 2 ) 1/e. Using Eqs. (|K|jl2|), and the l.h.s. 



of Eq. (prjf ), we see that the form factor k(M 2 ) is indeed finite through 0(a). It is also 
gauge invariant and independent of the Higgs mass, as the Higgs boson does not couple 
to the photon, and the Higgs vertices are suppressed by the Yukawa couplings of the 
leptons. Explicit expressions of k(M 2 ) can be found in ||,||. 

Because of the top loop contribution to A 7 ^(M|), however, the form factor does 
depend on the top mass, albeit very weaklyfj 

? / 9x a / 8 „\ , mi ^,M 2 , , 

fc( ^ = __,( 1 _. 4 .) ln ^ +0( _^ ) . ( 13 ) 

Numerically, the real part of the form factor is very close to one [[J. For m t = 175GeV, 
M w = 80.356GeV, and s 2 = 0.2316, 

sin 2 6 l e f f = s 2 (M 2 z ) + 4 x 1(T 5 (14) 

The unusual smallness of the one-loop electroweak radiative corrections (~ 2 x 10~ 4 ) 
is due to large cancellations between bosonic and fermionic contributions, which are 
individually over one order of magnitude larger 0. The imaginary part of fc(Mj), 
although comparatively large, gives a negligible contribution to the asymmetries 0. 

Eq. (|14D gives the one-loop prediction for sin 2 Q^ff in the SM: once s 2 is calculated 
from the inputs via the l.h.s. of Eqs. (|3]), sin 2 O^fl can be known with a great accuracy. 
2 For simplicity, s 2 is defined here without decoupling the top quark. That approach is considered 

loh. 
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Of course, s 2 depends very sensitively on M H and m t , through Af w and p, and has 
an uncertainty from the value of a(M z ) of about 0.1%. Moreover, unknown higher 
order contributions could add a significant uncertainty to its determination. In that 
respect, however, the situation for s 2 seems to be reasonably under control |IT|,p^] after 
the calculation of the three-loop QCD [[L3|] and of the two-loop 0(a 2 m 2 /M^) [|14| , |15 
corrections. 



Here I will concentrate on the analysis of higher order effects on Eq. (|14f ). Given 
the estimated intrinsic accuracy on s 2 , we can aim to a precision of about 10~ 4 on 
the form factor k(M 2 ). A first candidate are the QCD corrections to the hadronic 
loops. Fortunately, the scale of the process at hand is such that perturbative QCD can 
be reliably applied. As quark loops appear only in the amplitude of Fig.|l](b), QCD 
affects only the 7 — Z mixing as in the figure below, where the lowest order diagrams 
are displayed. Since the axial and the vector current do not mix even in presence of 
QCD, only the vector current correlator contributes to A^z(M 2 ), reducing the two-loop 
calculation to the old QED one [IB]. However, two- loop QCD corrections to the general 



electroweak vector boson and scalar self-energies are known exactly |T7j for any value of 



the quark masses and the external momentum, in terms of simple logs and dilogs. Even 
the three-loop massless quark contributions and the heavy top expansions of the same 
correlators are now available |13j, 18 . 



Unlike the case of the p parameter, however, QCD corrections are not important for 
the form factor k(M 2 ). Using |DJ, we see that the light quark loops, neglecting all mass 



effects, get a correction 

AMMl & ^A(^H j2 )^(|_ 4((3) + M ), (15) 

whose real part amounts to about —3 x 10~ 5 . Using the expansions in the third pa- 
per of Ref. |T7[, we can see that the top contribution of Eq. flT3] ) can be adjusted to 
include QCD corrections by the replacement log m 2 /M 2 — > log mf/M 2 (1 + a s {m t )/^) — 
(15/4) a 8 (mt)/ir. Numerically, these QCD corrections contribute about 7 x 1CT 5 to 
k(M 2 ). It seems that there is no reason to go beyond two-loop level in QCD. 

A relatively large two-loop effect which can be easily accounted for is given by the 
interference of the imaginary parts of A^z{M 2 ) and of the photon self-energy 
(Fig.0(b)). Indeed, light fermion loops induce large and additive imaginary parts in 
these two self-energies. We can resum the gauge-invariant fermionic part of the pho- 
ton self-energy multiplying A jZ (M 2 ) in Eq. (|10D by M 2 /(M 2 - A{fj(M 2 )\ m ) @, where 
]yfg is the fermionic part of the photon two-point function renormalized in MS. 
The interference between the two imaginary parts (the contribution of the real parts is 
negligible) increases the value of k(M 2 ) by 1.9xl0~ 4 , a quite large two-loop effect. The 
inclusion of this kind of contributions is necessary. 
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(6) 

Figure 2: Main two-loop electroweak topologies contributing to k. 

Another potential source of large higher order effects are heavy particles; it is well- 
known that they do not decouple in the SM, and that the corrections due to the heavy 
top are dominant in p at the one-loop level. The study of the leading effects of this origin 
is therefore very important, and is made possible by an expansion in powers of the heavy 
masses. Let us start briefly considering the possibility of heavy Higgs effects. 

At one loop, the radiative corrections to amplitudes without Higgs bosons among 
the external particles depend at most logarithmically on the Higgs mass, in the limit in 
which this mass is heavy. At two-loop the corrections can be at most quadratic in the 
Higgs mass. For the p parameter these quadratic terms were computed long ago, and 
more recently a complete analysis for electroweak observables has been performed [[19 



in 



the case where the Yukawa couplings are neglected. The result is that these corrections 
are completely negligible. Still, the Higgs boson could appear in the irreducible two- 
loop 7 — Z° mixing in association with its Yukawa coupling with the top. In this case, 
however, the amplitude is also proportional to m 2 , and is best considered in connection 
with a heavy top. 

Finally, we consider the two-loop effects of a heavy top quark. In the case of 
the amplitude at hand, we have seen that, at one-loop, the top contributes in a mild, 
logarithmic way (see Eq. (0)). At two loops, by analogy, we expect only 0(a 2 m 2 /M 2 7 ) 
effects, compared with the 0(a 2 m^/M^) appearing in Ap [pO] . Following the recent 
calculation of the 0(a 2 m 2 /M^) effects on Eq. (Jjty jTJJ], which turn out to be relatively 
important [HJ, the study of these quadratic two- loop contributions is necessary. The 
rest of this lecture is devoted to a detailed presentation of such calculation. 

At the two-loop level the diagrams contributing to the Z° decay amplitude and 
proportional to can be grouped into three categories: (i) products of one-loop Z° 
wave function renormalization by one-loop vertices and self-energies (Fig.^Jwith the first 
two wavy lines representing Z 0, s); (ii) products of one-loop 7 — Z° mixing amplitudes by 
vertices and self-energies (Fig|| with the second wavy line representing a photon); (iii) 
two-loop irreducible vertices and 7 — Z° self-energy, displayed in Figs.[l|(a,b). 

We are interested only in the two- loop contributions proportional to m 2 . It is easy 
to realize that only the last category meets this condition: we have to keep in mind 
that the one-loop vertices do not contain the top quark and that the 7 — Z° self-energy 
depends only logarithmically on mt (Eq. (]13D). We also notice that the wave function 
renormalization factor of the Z° contains only log mt in the heavy top limit. This can 
be understood on dimensional grounds, as this factor involves only the derivative of the 
Z° vacuum polarization function w.r.t. the external momentum. We therefore need to 
consider only the subset of two-loop irreducible vertices and 7 — Z° self-energies which 
can be proportional to m 2 . 

In addition to the two-loop diagrams, we need to include the m 2 part of the two- 
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loop MS counterterm of s 2 , and renormalize the parameters appearing in the one- loop 
correction Eq. flTOl) . With our choice of MS couplings and on-shell masses, when does 
renormalization of the one-loop amplitude introduce m 2 terms? From Eq. flHD we see 
that Ss 2 does not contain any such term, and the same is trivially true of 5a. As the 
one-loop amplitude is at most logarithmic in m t , we conclude that the renormalization 
of the couplings does not contribute at the level of our calculation. On the other hand, 
the mass renormalization of the top quark and of the vector bosons does contribute 
0(m 2 ) terms, because of the top dependence of the on-shell mass counterterms. The 
Goldstone bosons and the ghost fields that appear in the one-loop loops also need mass 
renormalization. In the t'Hooft-Feynman gauge employed here, this is implemented 
assigning them the corresponding vector boson on-shell mass counterterm. In the case 
of the Goldstone bosons, however, a spurious quartic m\ dependence must be cancelled 
by a tadpole contribution, as dictated by Ward identities. 

At 0(a 2 m 2 /M^), we can finally write the two- loop form factor as 



c A% +c - f -\M 2 z ) 



Ml 



_ T/(2+c.t) _ 
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5®i 



+ A 



extra i 



(16) 



where the term A extra is generated in Eq. (|T2| ) by reexpressing the ratio 2M^/M 2 m 
terms of c 2 , using the l.h.s. of Eq. (|3b|). As Ap is 0(am 2 /M^), this introduces a crucial 
term for the finiteness of our calculation. In n dimensions it is necessary to keep 0(e) 
terms in Ap, so that a finite part is generated: A extra = —N c g 4 / (167T 2 ) 2 m 2 /M 2 l/(2e) + 
finite part, where N c = 3 is the number of colors. The two-loop counterterm 5^s 2 is a 
crucial component of the calculation in |TJ|], can be checked with and reads 



5®i 



N c g 4 



m. 



1 /l 



'16tt 



2\2 



c 2 M 2 z 



13 



s' 



e V8 36 



(17) 



We now begin the study of the irreducible diagrams with the vertices. The presence of 
the top quark is only induced by loop insertions on the internal vector boson propagators. 
The relevant diagrams are shown in Fig||(a,b). Scalar loops are generally suppressed by 





Figure 3: Two-loop vertex diagrams containing the top. 

powers of the lepton mass, but at two-loop level the mixing + — W + takes place through 
the top in the diagram of Fig.^|(c). For what concerns the diagrams of Fig||(a,b), we 
observe that the leading mf contribution is cancelled by the mass renormalization of 
the corresponding vector boson. As the overall divergence of such graphs is logarithmic, 
and the quadratic subdivergence of the top loop is the same as in the counterterm 
5M 2 r = A WW {M^) fa ^4^,^(0), the use of on-shell masses in the one-loop vertices removes 



'■WW 

„2™2 



all 0(a m 2 /M^) vertex contribution at two-loop level. This shows a clear advantage of 
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the renormalization procedure adopted. The case of Fig.|3](c) is obviously very different 
and requires an explicit two-loop calculation, which I describe as an example of the heavy 
mass expansion method. 

A few kinematical simplifications allow us to write the amplitude in Fig.[|(c) as 



k 



p 



d n k d n p 

(27r) 2 "/i- 4e k 2 p 2 PlPK t PQ w 



18) 



where I used the abbreviations P™ 



P 



ML PK t 



[p - kf 



rn 



t ■ 



and PQ V 



(p — q) — Mw- First we notice that / is in fact a self-energy integral, because all the 
dependence on the external lepton momenta has dropped out. It can therefore be treated 
as a two-loop two-point function. 

In general, a two-loop two-point function cannot be expressed in terms of polylog- 
arithms or other known functions (of course in special cases, like for QCD corrections 
considered above, the two- loop integrals can be expressed in a compact way). The 
choice is then between numerical evaluation, which can be now very efficient P5fl , and 
approximate analytic evaluation. In particular, we are now interested in a heavy mass 
expansion p^] that enables us to expand the integral in Eq. ([18]) in powers of mt, keep- 
ing only the leading m\ term. The simple method here illustrated is suggestive of the 
strategy adopted in |TJ| . 

The integral / has a threshold at q 2 = M^, and is therefore not analytic at q 2 = M 2 . 
We cannot use a simple Taylor expansion in the external momentum in the approximation 
q 2 <C m 2 , but we have to resort to an asymptotic expansion. The identity 



(p - kf 



mi 



k 2 



+ 



2 p ■ k — p 2 



mi 



(k 2 



rn 2 ) 



[(p — k) 2 — m 



(19) 



helps us separate different regions of the domain of integration. Using Eq. ( |19D in /, we 
obtain a disconnected two-loop integral (product of two one- loop integrals), plus a two- 
loop integral with improved ultraviolet convergence in the k integration and improved 
infrared convergence in the p integration. Applying Eq. (|T9| ) twice to /, we obtain, up 
to disconnected integrals proportional to k ■ p that vanish by symmetry, 



d n k d n p 
(2vr) 2n /i- 4t 



2(k-pf 



k ■ p (2 k ■ p — p 



2\2 



k 2 p 2 K 2 P W PQ W k 2 K 2 p 2 P 2 PK t PQ, 



(20) 



with K t = k 2 — m 2 . The first integral can be easily reduced to the product of two well- 
known one-loop integrals, while the second is 0(l/m 2 ), as can be checked counting the 
degrees of infrared and ultraviolet divergence, and can be dropped. Finally, we are able 
to write the vertex contribution proportional to JUj as 



y(2+c.t.) 



Ncg^s 2 m 



(16tt 2 ) 2 MJ 



C0[M 2 z ,M w ,M w ] 



+ finite + O ( — 



(21) 



where C0[q 2 , mi, m<2\ is a convergent one-loop integral. The method followed in this ex- 
ample can be used in a similar way in the case of the irreducible 7 — Z° self-energies. 
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Schematically, we can summarize the effect of the heavy top expansion as in the figure 
above, where the one loop diagrams are evaluated exactly at the appropriate momentum 
q 2 = M 2 . The two-loop vacuum (q 2 = 0) integrals can always be reduced to a closed an- 
alytic expression containing dilogarithms following the strategy proposed, for ex., in p4 |. 



For what concerns the counterterm pieces, the procedure consists simply of expanding 
all masses in the one-loop A 1 z(M 2 ) according to m — > m — 5m, keeping all 0(e) terms. 
The combination of irreducible and counterterm pieces is 



A (2+c.t.) (M 2s _ N c9 smj 



2 r 



s 2 c 2 C0[M 2 z , M w , M w \ - \ + ^-s 2 ) - + finite 

8 36 / e 



(22) 



Putting together the various pieces of Eq. (|16|) the cancellation of divergences can be 
verified. In the notation of |L4]], the final result for Re k( 2 > top \M 2 ), in units A^ c (d/(47rs 2 )) 2 



m 2 /M^, reads 

Afra/onw,^ -211 + 24 M + 462 s 2 - 64 Ms 2 /3 s 2 \ nn ., rf , ,, „, 
Re AA;( 2 ^)(M 2 ) = — + f - - - J B0[Mj, M 2 W} M, 

(ht - 4) Vht (8 s 2 - 3) g(ht) (6 + 27 ht - 10 M 2 + M 3 ) (3 - 8 s 2 ) In ^ 



+ 108 108 (ht - 4) 

c\ ,<> /5 7^\, (M-l)(8s 2 -3) ^.M. 

nc 2 + — + — s 2 )\nzt + - . . — -#(—)■ (23) 

6 U 18 / 18 (4 -fa) fa V 4 ; v ; 

Numerically, this correction is negligible. It grows with M H , reaching about 2.5 xlO -5 
for M H = 700GeV, well below our precision goal. However, it must be emphasized that 
this result is strictly related to the renormalized parameters that have been chosen here, 
and to the way the one-loop result has been expressed. When all the two-loop effects 
considered here are included, and using the same inputs, Eq. fll4|) is modified into 

sin 2 9 l e f f = s 2 (M 2 z ) + 1.0 x 10~ 4 . (24) 

In summary, I have reviewed the radiative corrections entering the relation between 
the effective leptonic sine and the MS parameter sm 2 9 w (M 2 ). The main higher order 
effects have been considered, and a new calculation of the 0(ct 2 m 2 /M^) contributions has 
been illustrated in detail. All these corrections are very small and we can conclude that 
in the case at hand higher order effects are unlikely to be a major source of uncertainty. 

I am grateful to Marek Jezabek and to the organizers of the School for the kind 
invitation, the hospitality, and the excellent organization. The results of the two-loop 
heavy top calculation have been obtained in collaboration with G. Degrassi. 
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